Astronomy & 
July 3, 2008 



Astrophysics manuscript no. abridged 



© ESO 2008 



NGC300 X-l and IC10 X-l: a new breed of black hole binary? 



R. Barnard, J. S. Clark, U. C. Kolb 



00 

O 

o 

(N 



S 1 
6 

CO 



> 

O 

q 
o 

oo 

o 



% 



The Department of Physics and Astronomy, The Open University, Walton Hall, Milton Keynes, MK7 6BT, U.K. 

Received / Accepted 

ABSTRACT 

Context. [ABSTRACT ABRIDGED] IC10 X-l has recently been confirmed as a black hole (BH) + Wolf-Rayet (WR) 

X-ray binary, and NGC300 X-l is thought to be. IC10 X-l and NGC300 X-l have similar X-ray properties, with 

luminosities ^lO* 3 erg s _1 , and orbital periods ~30 hr. 

Aims. We investigate similarities between these systems, as well as differences between them and the known Galactic 

BH binaries. 

Methods. We have examined XMM-Newton observations of NGC300 X-l and IC10 X-l. For each observation, we 

extracted lightcurves and spectra; power density spectra (PDS) were constructed from the lightcurves, and the X-ray 

emission spectra were modeled. 

Results. Each source exhibits PDS that are characteristic of turbulence in wind accretion or high state X-ray binaries 

(XBs). In this state, Galactic XBs with known BH primaries have soft, thermal emission; however the emission spectra 

of NGC300 X-l and IC10 X-l in the XMM-Newton observations are predominantly non-thermal. The remarkable 

similarity between the behaviour of NGC300 X-l in Observation 1 and that of IC10 X-l lends strong evidence for 

NGC300 X-l being a BH+WR binary. 

Conclusions. The unusual spectra of NGC300 X-l and IC10 X-l may be due to these systems existing in a persistently 

high state, whereas all known BH LMXBs are transient. BH XBs in a persistent high state could retain their corona, 

and hence exhibit a large non-thermal component. LMC X-l is a BH XB that has only been observed in the high 

state, and its spectrum is remarkably similar to those of our targets. We therefore classify NGC300 X-l, IC10 X-l and 

perhaps LMC X-l as a new breed of BH XB, defined by their persistently high accretion rates and consequent stable 

disc configuration and corona. This scenario may also explain the lack of ultraluminous X-ray sources in the canonical 

soft state. 

Key words. X-rays: binaries - X-rays: general - Galaxies: individual: NGC 300 - Galaxies: individual: IC 10 



1. Introduction 

Characterised by high temperatures and large mass loss 
rates, Wolf Rayet (WR) stars represent the final stage of 
stellar evolution prior to supernova for massive stars (>25 
M Q ; ICrowtherl |2007l ). Their spectra are characterised by 
strong, broad emission lines, which reveal their chemically 
evolved nature due to, successively, the CNO (WN stars) 
and triple-a (WC stars) processes. These products of nu- 
clear burning become visible as the hydrogen rich mantle 
of the O star progenitor is stripped away. For single stars 
this is accomplished via the stellar wind of the progenitor; 
however binary interaction may also contribute to this pro- 
cess. Indeed, if a stellar merger may be avoided, interaction 
between the two components of a high mass X-ray binary 
(OB star and compact object) is thought to lead to a short 
period WR+compact object binary (WR+co). 

Empirically such systems are exceedingly rare; to 
date we know of only three possible WR+co sys tems: 
Cygnus X-3 in our Galaxy (iLqmmen et all |2005|), IC 
10 X-l dClark fc Crowtherl I2004t iBauer fc Brandt! 120041) 



and most recently, NGC 300 X-l (jCarpano et all 120071: 
Crowther et alj|2007h . Population synthesis carried out by 
Lommen et al.l (|200 5) suggests that ~1 WR + BH binary 
and ~1 WR + neutron star (NS) binary, may exist as bright 
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X-ray sources in a galaxy like our own; this is certainly in 
keeping with the rarity of candidates. 

iPrestwich et al.l (|2007f ) re-analysed archival Gemini 
Multi-Object Spectrograph spectra of the WR star as- 
sociated with IC10 X-l, and found evidence for a black 
hole (BH) primary, with a plausible BH mass range of 
23-34 Mfry This result was subsequently confirmed by 
ISilverman fc Filippenkol (|2008T ). who obtained spectra over 
five nights of observations with the Keck-I 10m telescope. 
They report a BH mass of 20-35 M for the compact ob- 
ject, making it the most massive known stellar BH. NGC300 
X-l has not yet been optically identified as a WR+BH 
binary; however, in this work we find striking similari- 
ties between the X-ray properties of NGC300 X-l and the 
confirmed WR+BH binary IC10 X-l. We will show that 
NGC300 X-l is likely to be a WR+BH binary also. 

WR+co binaries are just one subset of high mass X-ray 
binaries (HMXBs), and it is possible that their temporal 
and spectral characteristics closely resemble other types of 
HMXBs. Hence it is instructive to briefly review the perti- 
nent properties of HMXBs in general. HMXBs are well clas- 
sified by the nature of the donor star (see e.g. IWhite et al.l 
1995], and references within). 

Around 60% of Galactic HMXBs include a Be donor 
()Liu et al.l |2007| ); these have elliptical orbits lasting tens 
to hundreds of days, only accreting near periastron. As 
a result, Be HMXBs are transient sources, generally pow- 
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ered by Bondi-Hoyle accretion (Bondi & Hoyle 1944) from 
a dec retion disc around the Be star (jOkazaki fc Negueruelal 
120011) . Several Bondi-Hoyle accreting HMXBs with a NS 
accretor exhibit rapid aperiodic variability in X-ray inten- 
sity with amplitude ~10-30%; their power density spectra 
(PDS) resemble a power law with spectral index 7 ~1 ( 
powe r oc z/~ 7 where v is the frequency, see e.g. iNagasd 
119891 iBelloni fc Hasingerlll99d Ivan der Klid 11995ft . These 
works suggested that such variability could be present in 
all accretion-powered pulsars. However, an inspection of 
the cat alogues of HMXBs in the Galaxy and Magellanic 
clouds (jLiuet al.l [20051120071 ) revealed that only 16 out of 
242 HMXBs have published records of aperiodic variabil- 
ity. HMXBs that have exhibited the variability described 
above have power law emission spectra with spectral in- 
dex, T, ~0.5-1.8. There arc no known BH HMXBs that 
accrete in this way. 

A further ~32% of Galac tic HMXBs contain supergiant 
(SG) donors (|Liu et al.ll2007l ); the orbital periods range over 
1.4-41 days, and the donor stars are (or are close to) filling 
their Roche lobe (see lKaper et al.ll2004 for a review). The 
compact object in SG HMXBs is continuously accreting 
from the donor's wind and can be either wind-fed via Bondi- 
Hoyle accretion (with luminositi es ~10 35 -10 36 erg s -1 ) or 
disc-fed, reaching ~10 38 erg s _1 ()Kaper et al,ll2004ft . 

Studies of Galactic disc-accreting XBs have shown that 
their spectral and variability properties depend as much on 
the accretion rate as the nature of the accretor (NS or BH). 
At low accretion rates, the emission spectra and temporal 
variability of disc-accreting binaries are remarkably similar 
(|van der K lis 1994). However, at certain, higher, accretion 
rates there are characteristic differences that may allow one 
to determine the nature of the primary; the PDS are well 
described by a power law with index 7 ~1, but the emis- 
sion spectra for NS and BH XBs are very different. For a 
disc- accreting XB with a NS primary, the emission corre- 
sponding to the power law PDS is non-thermal, while disc- 
accreting BH XBs exhibiting the same PDS hav e a ther- 
mal X-ray emission spect rum (jvan der Klisl ll994). Indeed, 
iDone fc Gierlinskil (|2004ft examined - 1 Tbyte of RXTE 
archival data on Galactic XBs, and found a region in colour- 
colour space that is uniquely associated with BH XBs in 
this high, soft state. Hence, we may learn a great deal 
about disc-accreting binaries from their X-ray variability 
and emission spectra. 

IC10 
X-ray 
200 
200- 
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Table 1. Journal of XMM-Newton observations of 
NGC300 X-l and IC10 X-l. For each observation, we pro- 
vide the target, date, revolution number, exposure and good 
time interval. 



Object 


Date 


Rev 


Exp 


GT 


NGC300 X-l 


2000 Dec 26 


0192 


32 ks 


25 ks 


NGC300 X-l 


2001 Jan 01 


0195 


40 ks 


40 ks 


NGC300 X-l 


2005 May 22 


0998 


35 ks 


25 ks 


NGC300 X-l 


2005 Nov 25 


1092 


35 ks 


35 ks 


IC10 X-l 


2003 Jul 03 


0653 


44 ks 


28 ks 



IC10 X-l, in order to compare their properties. In partic- 
ular, we examined the PDS from the X-ray lightcurves of 
these systems for the first time. In Sect. [2] we describe the 
observations and data analysis, then provide the results 
from NGC 300 X-l and IC10 X-l in turn in Sect. [3J We 
discuss our findings in Sect. [H and consider four scenarios 
for NGC300 X-l and IC10 X-l: wind accretion onto a 
neutron star, disc accretion onto a neutron star, wind 
accretion onto a black hole and disc accretion onto a black 
hole. We then compare the properties of these sources 
with those of the ultralurninous X-ray sources (ULXs), 
many of which are thought to be HMXBs. Finally we draw 
our conclusions in Sect. El 



2. Observations and data analysis 

Four XMM-Newton observations have been made of NGC 
300 X-l, and one of IC10 X-l; a journal of observations is 
provided in Table [TJ For our analysis we used the XMM- 
Newton SAS version 7.0, and the FTOOLS suite, version 
5.3.1. For each observation, we filtered out intervals of high 
background (flaring), using the criteria recommended by 
the SAS team. We note that no flaring occurred during the 
2001, January or 2005, November observations of NGC300 
X-l. 

We then extracted pn and MOS 0.3-10 keV lightcurves 
and spectra from a circular region centred on the source, 
along with corresponding response files. Background re- 
gions were then chosen, and lightcurves and spectra were 
obtained from these regions for pn and MOS also. 

The pn and MOS source and background lightcurves 
were co-added, after careful synchronisation; XMM 
lightcurves are non-synchronised by default, leading to ar- 
tificia l variability if not treated properly ([Barnard et al.l 
120071 ). Background-subtracted lightcurves were analysed for 
variability, and PDS were made from the combined EPIC 
source lightcurves (background not subtracted). These 
lightcurves were averaged over several intervals of 1024 bins, 
with 5.2 s binning and geometric grouping; the sampled fre- 
quency range was ~0. 0002-0.1 Hz. 

The source and background emission spectra from the 
two MOS cameras were combined with the FTOOL math- 
pha to give MOS1+MOS2 source and background spectra, 
and the corresponding response matrices and ancillary re- 
sponse files were combined also, with addrmf and addarf. 
We modeled the pn and combined MOS spectra simultane- 
ously using XSPEC 11.3, with a constant of normalisation 
to account for differences in the pn and MOS responses. 



3. Results 

We discuss our analysis of the variability and X-ray spectra 
from NGC300 X-l and IC 10 X-l in the following sections. 



3.1. Variability of NGC300 X-l and IC10 X-l 

The lightcurves from the four XMM -Newton observations 
of NGC300 X-l are well described bv lCaroano eTaD lj2007ft 
and will not be discussed further here. Each observation re- 
veals strong variability, with PDS that are well described 
by power laws with 7 ~1. We present the 2005, November 
PDS in Fig. [TJ The axes are log scaled, and normalised to 
give r.m.s. 2 variability; the expected noise is subtracted. We 
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Fig. 1. PDS of the combined pn+MOS 0.3-10 keV 
lightcurve of NGC300 X-l from the 2005 November XMM- 
Newton observation (Obs. 4). The axes are log scaled, and 
the y axis is normalised to give r.m.s. 2 variability. The ex- 
pected white noise level (7.48) is subtracted. A power law 
fit to the PDS is shown, with 7 = 1; the best fit x 2 /d.o.f 
= 26/24, with a good fit probability of 0.35. Such a PDS is 
characteristic of a disc-accreting X-ray binary with a high 
accretion rate. 
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Fig. 2. Combined EPIC lightcurve of IC10 X-l from the 
2003 XMM-Newton observation (top); the binning is 100 
s. The lower panel shows the ratio of 2.5-10 keV photons 
to 0.3-2.5 keV photons (hardness ratio) for each time bin. 
The lightcurve is clearly variable, and there is evidence for 
variability in the hardness ratio. 



note that no background flaring occurred in this observa- 
tion, so there is no question of the PDS being an artefact of 
flaring or background filtering. This variability is certainly 
significant, as fitting the PDS with zero power yields a \ 2 
of 74 for 24 degrees of freedom (dof). We present the r.m.s. 
variability and best fit 7 for each observation of NGC300 
X-l in Tabled We note that the PDS from Obs. 3 is dif- 
ferent to the PDS from the other observations; this differ- 
ence is likely due to the eclipse that occurs in Obs. 3 (see 
ICarpano et ai]|2007t) . 



We present the 0.3-10 keV, combined pn+MOS 
lightcurve of IC10 X-l in the top panel of Fig. [21 the back- 
ground has been subtracted, and intervals of high back- 
ground have been removed. The lightcurve has 100 s reso- 
lution. The system is highly variable throughout the obser- 
vation, with a large intensity dip near the beginning. This 
dip may be intrinsic to the X-ray source, or due to an in- 
crease in line-of-sight absorption, such as an eclipse. The 
bottom panel shows the hardness ratio for each time bin, 
defined as the ratio of 2.5-10 keV counts to 0.3-2.5 keV 
counts. The hardness ratio is variable at the 3tr level, and 
appears to be higher during the intensity dip at the start of 
the observation, hinting at photo-electric absorption; how- 
ever, we cannot determine the phase of this intensity dip 
with the current orbital ephemeris. The variability is also 
consistent with stochastic variations in a source with a PDS 
described by 7 ~1. 

The PDS created from this lightcurve is presented in 
Fig.[3J as before, the y-axis is normalised to show the r.m.s. 2 
power, and the expected noise is subtracted. This PDS is 
also acceptably fitted by a power law with 7 = 1, with 
x 2 /dof = 35/24. As with NGC300 X-l, the r.m.s. variability 
is greater than for Galactic LMXBs, at 18±2%, excluding 
the large intensity dip. 

The PDS observed from NGC300 X-l and IC10 X-l 
are characteristic of disc-fed XBs at high accretion rates, 
albeit with fractional r.m.s. variabilities ~5— 10 times highe r 
than observed in Galactic XBs (see e.g. Ivan der Klis! [l995'). 
However, they are also similar to the variability observed 
in certain wind-fed HMXBs. Hence, we must examine the 
emission spectra of NGC300 X-l and IC10 X-l if we are to 
determine their natures. 

3.2. Modeling the spectra of NGC300 X-l and IC10 X-l 
3.2.1. NGC300 X-l 



ICarpano et al.l (|2007h modeled the spectra from the four 
observations of NGC300 X-l together. They used an ab- 
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Fig. 3. Combined EPIC PDS of IC10 X-l from the 2003 
XMM-Newton observation. The axes are log-scaled, and the 
y-axis is normalised to give r.m.s. 2 variability. The expected 
noise level (3.49) is subtracted. The best fit power law with 
7 = 1 is shown, for % 2 /dof = 35/24. 
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1 2 

Energy (keV) 

Fig. 4. Top panel: Best fit model to the Obs. 1 pn (black) 
and combined MOS (grey) 0.3-10 keV spectra of NGC300 
X-l using a power law with T = 2.44 and a 0.95 keV 
Gaussian line with width 0.07 keV; x 2 /dof = 100/60, which 
is unacceptable Bottom panel: x 2 residuals, showing a hard 
excess; this shows that a flatter power law and/or a second 
emission component is required. The spectrum of NGC300 
X-l from Obs. 1 is certa i nly di fferent from the model ob- 
tained by ICarpano et all (|2007t ) when fitting spectra from 
all four observations simultaneously. 



1 2 

Energy (keV) 

Fig. 5. Simultaneous fits to the pn (black) and combined 
MOS (grey) spectra of IC10 X-l from the 2003 XMM- 
Newton observation with the best fit emission model; the 
first 20 ks of the observation was excluded. The top panel 
shows the spectra modeled with a power law with photon 
index 3.5±0.2 and a 0.86±0.03 keV keV blackbody, suffer- 
ing absorption equivalent to 8.6±0.4xl0 21 H atom cm~ 2 , 
with x 2 — 849 for 840 degrees of freedom (good fit proba- 
bility 0.4). The bottom panel shows the x 2 residuals. 



sorbed, two-component model, consisting of a power law 
with photon index T, i.e. N(E) oc E~ T , where N is the 
number of photons at energy E, and an emission line. They 
imposed identical values for T, line energy and line width 
for the four observations, with these values free to vary; the 
absorption and normalizations for each component were in- 
dividually varied for each observation. For their best model, 
T = 2.44±0.03 , the line energy was 0.94±0.01 keV, with a 
width of 0.07±0.01 keV. They quote a reduced x 2 of 1.15 
for 1033 dof; however, this equates to x 2 /dof = 1188/1033, 
which is rejected at a 99.95% level and is therefore unac- 
ceptable. We therefore examined the spectra of NGC300 
X-l more carefully. 

For each observation of NGC300 X-l, we modeled the 
pn and combined MOS spectra simultaneously, with each 
model consisting of a power law + emission line, suffer- 
ing absorption and including a constant of normalization 
to account for differences in pn and MOS calibration. The 
resulting best fits are shown in Table [3J We found the spec- 
tr a from Obs. 2-4 to be consistent with the model reported 
bv lCarpano et al.l (|2007l ). However, the best fit to the Obs. 
1 spectra yielded x 2 /dof = 100/60, an unacceptable fit. We 
show this best fit in Fig. [4] along with x 2 residuals. We 
see that the model systematically underestimates the emis- 
sion at higher energies, suggestive of an additional emission 
component. Hence, we fitted all NGC300 X-l spectra with 
absorbed blackbody + power law models and also absorbed 
disc blackbody + power law models, commonly observed in 
NS and BH XBs respectively. The best fits for these models 
are also provided in Table [3J 0.3-10 keV lu minosities are 
provi ded, assuming a distance of 1.88 Mpc (|Gieren et al.l 
2005). We see that both two-component models provided 
good fits to the Obs. 1 and Obs. 3 spectra. The Obs. 3 
spectrum is well described by all three models; however, the 
resulting luminosities differ by factor of two (~5er); hence 
the nature of the Obs. 3. emission spectrum is uncertain. 



Neither the Obs. 2 nor Obs. 4 spectra were successfully de- 
scribed by any thermal + power law model; hence NGC3000 
X-l exhibits at least two spectral states: one exhibited in 
Obs. 1 and the other exhibited in Obs. 2 and 4. 

Finally, we fitted the spectra with the hard power 
law emission that is characteristic of Bondi-Hoyle accretion 
onto a neutron star. Our literature review shows that the 
wind-accreting HMXBs with PDS similar to those observed 
from NGC300 X-l all have spectra that are well described 
by a power law with T <1.8. Hence, we modeled the spectra 
from each observation with power law models, restricting 
T to be <1.8. We found the best fit x 2 /dof to be ^5 in 
each case. Hence NGC300 X-l is unlikely to be powered by 
Bondi-Hoyle accretion onto a neutron star. 

3.2.2. IC10 X-l 

We first attempted to simultaneously fit the pn and com- 
bined MOS spectra of IC10 X-l with various emission mod- 
els, but failed to obtain a good fit. This is likely to be due 
to the fact that IC10 X-l varies in colour and intensity 



Table 2. Fits to the 0.3-10 keV PDS of NGC300 X-l 
from Obs. 1-4. We give the observation number, as well 
the best fit spectral index with the x 2 /dof. We also give 
the x 2 /dof for zero power. Finally we give the fractional 
r.m.s. variability. We note that the PDS for Obs. 3 is more 
complex than for the other observations, due to an eclipse. 



Obs. 



X 2 /dot x 2 / dof (0) Frac. r.m.s. (%) 



Obs. 1 


1.30(12) 


17/23 


57/25 


50(3) 


Obs. 2 


1.3(2) 


27/23 


66/25 


44(2) 


Obs. 3 


1.8(2) 


56/23 


91/25 


38(2) 


Obs. 4 


1.05(11) 


23/23 


74/23 


22(2) 
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Table 3. Fits to the spectra of NGC300 X-l from Obs. 1-4, and of IC10 X-l. The first model is an absorbed power 
law plus Gaussian emission line, as favoured by Carpano et al. (2007a). For this model, we give absorption, normalised 
to 10 22 H atom cm -2 , spectral index, line energy and width (in keV), and x 2 /dof, with good fit probability in square 
brackets. We only apply the Gaussian to NGC300 X-l. We also give the 0.3-10 keV luminosity normalised to 10 37 erg 
s _1 . We next provide the best fits for a disk blackbody + power law model that is often used when describing the high 
state of disc-accreting BH XBs. We give the absorption, temperature of the disk blackbody, spectral index of the power 
law component and % 2 . We then give 0.3-10 keV luminosity as before, and the fractional contribution of the power law 
model. 









NGC300 X-l 






IC10 X-l 


Model 


Parameter 


Obs. 1 


Obs. 2 


Obs. 3 


Obs. 4 




PL + Line 


"ii 


0.036 


0.067(12) 


0.052(9) 


0.07(2) 


0.75(3) 




r 


2.46(6) 


2.44(7) 


2.44 


2.53(8) 


2.56(6) 




-BLinc 


0.95 


0.93(2) 


0.95(2) 


0.95(2) 






CT 


0.07 


0.08(2) 


0.07 


0.06(3) 






X 2 /dof 


100/60 [9E-4] 


149/154 [0.60] 


50/60 [0.82] 


117/145 [0.05] 


794/652 [1E-4] 




-^0.3-lOkeV 


11.5(13) 


28.4(12) 


20(2) 


34.8(16) 


39(3) 


BB+PL 


"II 


0.11(4) 


0.130(12) 


0.22(5) 


0.17(3) 


0.92(10) 




kT 


1.02(14) 


1.02(17) 


0.92(12) 


1.08(15) 


0.82(3) 




r 


3.41(6) 


3.01(16) 


3.8(3) 


3.26(6) 


3.9(5) 




X 2 /dof 


52/58 [0.67] 


211/155 [2E-3] 


59/58(0.43] 


221/148(1. 0E-4] 


697/650 [0.09] 




^0.3-lOkeV 


16.9(17) 


39(2) 


49(6) 


59(3) 


91(19) 




Frac. PL 


0.80(11) 


0.93(14) 


0.91(18) 


0.92(6) 


0.93(3) 


DISKBB+PL 


UH 


0.12(4) 


0.13(4) 


0.22(9) 


0.17(4) 


0.84(4) 




klin 


1.8(4) 


1.9(6) 


1.5(4) 


2.0(4) 


1.19(6) 




r 


3.6(7) 


3.1(4) 


4.1(7) 


3.4(3) 


3.9 




X 2 /dof 


56/58 [0.54] 


213/155 [1.3E-3] 


62/58(0.33] 


226/148[3.0E-5] 


710/651 [0.05] 




^0.3-lOkeV 


19(3) 


33(25) 


49(6) 


53(5) 


72(9) 




Frac. PL 


0.89(17) 


0.6(7) 


0.91(18) 


0.93(11) 


0.86(19) 



during the observation, mixing the spectral states so that 
a physical fit is impossible. 

In light of this, we excluded the first 20 ks of the obser- 
vation, where the intensity and colour of IC10 X-l is unsta- 
ble. We modeled the spectrum of the remaining ~22 ks with 
power law, blackbody + power law, and disc blackbody + 
power law emission models, all suffering absorption. The 
best fits are presented in Table [3l quoted lumin osities as- 
sume a distance of 740 kpc (jDemers et al.l 120041 ). The disc 
blackbody + power law model favoured an unreasonably 
high r, although T was not well constrained. Hence, we 
fixed it to 3.9, the best fit T for the blackbody + power law 
model. No simple power law model was able to fit, but both 
two-component models provided acceptable fits. However, 
since we know that IC10 X-l contains a BH, we favour the 
disc blackbody + power law model. We present the best fit 
to the IC10 X-l spectrum with the disc blackbody + power 
law model in Fig. [5] 

IWang et alj (|2005l ) modeled pn and MOS spectra 
from the entire XMM-Newton observation, mixing spectral 
states and thereby producing a spectrum that is not repre- 
sentative of either state. As a result, they required a metal 
abundance <0.01 Solar to achieve good fits; they remark 
that such a low metallicity is probably not physical, but 
attribute it to poor calibration at low energies rather than 
to spectral mixing. Their preferred model is a self-consistent 
Comptonised multi-colour disc blackbody where an accre- 
tion disc emits a multi-temperature disc blackbody, which 
is Comptonised by a spherically s ymmetric co r ona o f hot 
electrons around the disc. Hence IWang et al.l (|2005f ) also 
favour a non-thermal emission model for IC10 X-l. They 
estimated the mass of the black hole by assuming that the 



inner disc radius corresponds to the innermost stable orbit, 
after correcting for relativistic effects. They obtained a BH 
mass range of ^4 M Q for a stationary BH, up to ^30 M© 
for a BH with extreme spin; this range overla ps the mass 
range found bv lSilverman fc Filippenka (|2008f h 

4. Discussion 

We have examined the variability and emission spectra of 
one BH+WR candidate (NGC300 X-l), and the newly con- 
firmed BH+WR system IC10 X-l. The 0.3-10 keV spec- 
trum of Obs. 1 of NGC300 X- 1 is strikingly similar to that of 
IC10 X-l. This lends strong support to the case for NGC300 
X-l being a BH+WR system. Both systems exhibit vari- 
ability that scales with frequency, v, as v , identified with 
BH XBs in the high state. However, their 0.3-10 keV emis- 
sion spectra appear to be ~90% non-thermal, while the 
canonical high state BH spectrum is thermally dominated. 
We therefore suggest that these systems exist in a previ- 
ously unrecognised state. 



4.1. Accretion scenarios 

We next discuss different accretion scenarios. IC10 X-l has 
been confirmed as a WR+BH binary, hence we first discuss 
wind and disc accretion onto a BH. NGC300 X-l has an un- 
known accretor, and may contain a neutron star; therefore, 
we also discuss wind and disc accretion onto a NS. 

Since the emission from Bondi-Hoyle accreting neutron 
stars appears to originate on or near the surface of the 
neutron star, it is unclear what X-ray emission could be 
expected from Bondi-Hoyle accretion onto a black hole. 
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Indeed, there are no known black hole Be XBs, and all 
known black hole supergiant XBs exhibit accretion discs. 
Furthermore, the modeling of Bondi-Hoyle accretion onto 
black holes requires relativistic, thre e-dimensional m agne- 
tohydrodynamical modeling (see. e.g. I Font et alJll999 . and 
references within) . These models have not yet produced un- 
ambiguous predictions for the variability or emission mech- 
anisms expected from such systems. Hence, there are no ob- 
servational or theoretical constraints on the spectral shape 
or variability exhibited by HMXBs powered by Bondi-Hoyle 
accretion onto a BH. 

The short orbital periods (~30 hr) and high X-ray lu- 
minosities of NGC300 X-l and IC10 X-l rese mble those of 
disc-accreting SG HMXBs (|Kaper et al.ll2004D; we therefore 



consid er them likely disc-accretors also. iPrestwich et al.l 
(|2007f ) calculated the range of orbital periods that permit 
Roche lobe overflow for IC10 X-l. They obtained periods 
of ^2.5-3 hr, and rejected disc accretion for IC10 X-l (this 
also applies to NGC300 X-l, as it has a similar orbital pe- 
riod to IC10 X-l). However, it is entirely possible for the 
winds of the Wolf-Rayet stars to power disc accretion. The 
three known black hole HMXBs Cygnus X-l, LMC X-l and 
LMC X-3 are all disc-accreting, with orb ital periods of 5.60 
d, 4.22 d and 1.70 d respectively (see e.g. lLewin et al.|[l995l 
and references within). Of these, only LMC X-3 is thought 
to be Roche lobe filling. NGC300 X-l and IC10 X-l have 
shorter orbital periods than any of these systems, so disc 
accretion onto a BH in these systems is entirely plausible. 
As discussed in Sect. 3, our spectral modeling of 
NGC300 X-l allows us to rule out Bondi-Hoyle (wind) ac- 
cretion onto a NS, as power law spectral fits with T < 1.8 
yield x 2 /&oi >5 for all observations. However, the emis- 
sion and variability of NGC300 X-l are in keeping with a 
disc-fed NS XB. The donor star may be the WR star, in 
which case, the disc would likely be wind- fed. Alternatively, 
NGC300 X-l could be a bright NS LMXB which is merely 
coincident with the WR. 



4.2. A new black hole state? 

We prefer disc accretion scenarios for IC10 X-l and 
NGC300 X-l, although there is no independent evidence for 
a disc, such as radio jets. All known BH XBs are disc accret- 
ing, yet the properties exhibited by IC 10 X-l and NGC300 
X-l are as yet unrecognised. Their observed emission spec- 
tra most clos ely resemble the ve r y high (steep power law, 
SPL) state (jvan der Klisl 119951 : IMcClintock fe Remillardl 
2006). However, the PDS of sources in the SPL state 
are characterised by broken power laws, where 7 changes 
from ^0 at low frequencies to ~1 at high frequen- 
cies; often quasi-periodic oscilla tions (QPOs) are observed 
(|McClintock fc Remillardll2006h . 

We note that Cyg X-3, the other WR+co candidate, ex- 
hibits similar timing behaviour to NGC300 X-l and IC10 
X-l. Some believe that scattering in the wind suppresse s 
high frequency variability (see e.g. iKitamoto et al. 1992), 
possibly masking a SPL PDS. Hence, one might wonder if 
NGC300 X-l (in Observation 1) or IC10 X-l are disguising 
SPL states. However, the r.m.s. variability of NGC300 X-l 
was already 44% in Observation 1, while the r.m.s. variabil- 
ity of IC10 X-l was 18%. Hence they are already consider- 
ably more variable than BH LMXBs in the very high state 
(|van der Klislll995l ). making this hypothesis unlikely. 



Instead, we believe that the observed behaviour maybe 
due to the status of the corona. In NS LMXBs, the 
self-irradiation of the disc is more efficient than in BH 
LMXBs, hence NS LMXBs have hot, stable discs , while 
all BH LMXBs are transient (jKing et all Il997t) ; they 
spend the majority of their time in a quiescent state, 
with outbursts generally last ing a few months (see e.g. 
IMcClintock fc Rem illard 2006) . During the rise of the out- 
burst, a BH LMXB will go from a low accretion rate, hard 
state with associated radio j ets, to a high accretion rate, 
soft state with no jets (see e.g. lHoman fe Bellonil l2005. and 
references within) . It is therefore possible that the source of 
the hard, non-thermal emission, and also of the jets in BH 
LMXB s (e.g. the c o rona) is ejected during this transition; 
indeed iGallo et al.l (|2004h witness a large radio flare that 
coincided with a hard to soft transition in GX 339—4. 

By contrast, the three known Galactic BH HMXBs are 
all persistently bright. The HMXB accretion discs are ex- 
pected to be small, and hence will be more easily kept hot 
at the outer edge, making them stable. 

Therefore, the behaviour of IC10 X-l and NGC300 X-l 
may be explained by stable disc accretion. A small fraction 
of the ~ 10~ 5 — 10~ 4 M Q yr _1 wind from a WR companion 
leaving the vicinity of the LI point could provide enough 
mass, with sufficient angular momentum, to sustain an ac- 
cretion disc that is persistently in the high state. In this 
case, the corona could be retained, and the emission would 
then be dominated by a non-thermal component due to 
inverse-Comptonisation of cool photons on hot electrons, 
as seen in NS XBs at high accretion rate. 

Observations of the Galactic BH HMXBs support this 
hypothesis. LMC X-l has never been observed in th e 
low/hard stat e (see e.g. IWilms et al.ll2001l:lYao et al.ll2005l ). 
iNowak et al.l (|2001h report "high state" (i.e. power oc 
y- 1 ) PDS from R XTE observations of LMC X-l, while 
I Wilms et al.l ([200 1[ ) find its emission spectrum to be domi- 
nated by a power law with photon index <;3, with a 0.8-1.1 
keV blackbody component. This emission spectrum is re- 
markably similar to our best fit models for NGC300 X-l in 
Obs. 1 , and to IC10 X-l. In contrast, Cyg X-l and LMC X- 
3 exhibit transitions between the low/hard and high states, 
and exhibit canonical spectra (|McClintock fc Remillardl 
2006, and references within). 

4.3. Comparison with ULXs 

It is interesting to compare these systems with the ULXs, 
as many ULXs are thought to be HMXBs. In their re- 
view of black hole m asses and spectral states in ULXs, 
ISoria fc Kuncid ((2(308) find that only ~10% of sources with 
luminosities above 10 erg s _1 exhibit high/soft (i.e. ther- 
mally dominated) spectra, with almost all ULXs instead 
dominated by a broad "power law-like" component that 
contributes ~90% of the flux. Such characteristics are re- 
markably similar to those of NGC300 X-l and IC10 X-l, 
where the power law contributes ~80-90% of the flux. It 
is therefore possible that these ULXs could also be pere- 
sistently bright with stable coronae. It is therefore instruc- 
tive to study the spectral and, crucially, timing properties 
of ULXs. 

IWinter et al.l (|2006l ) conducted a survey of all point 
sources >10 38 erg s _1 in XMM-Newton observations of 32 
nearby galaxies. They modeled the spectra of ^100 X-ray 
sources, each with >400 counts. Ten sources exhibited spec- 
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tra consistent with NGC300 X-l in Obs. 2 and 4, including 
three ULXs. A further 23 sources exhibited spectra consis- 
tent with Obs. 1 of NGC300 X-l and IC10 X-l, including 
10 ULXs. Hence the processes in NGC300 X-l and IC10 
X-l could also occur in some ULXs. 

While most of the known ULXs are too faint to al- 
low fruitful analysis of the PDS, several of the nearest 
ULXs are as bright as, or brighter than, NGC300 X-l 
and IC10 X-l. A full analysis of the timing properties 
of nearby ULXs is beyond the scope of this work, and 
is the subject of a forthcoming paper. However, we con- 
ducted a brief literature review of ULX variability. Several 
ULX s exhibit Q POs, e.g. NGC5408 X-l (IStrohmaver et all 
I2007IL M82X-1 (IStrohmaver fe Mushotzkvll2003l) and HoLX 
X-l (jDewangan et al.ll2006D . IStrohmaver et all (|2007l ) the 
NGC5408 X-l PDS to be characterised by a broken power 
law; this PDS, along with the emission spectrum (domi- 
nated by a power law with F ^2.6) is reminiscent of the 
SPL state of Galactic black holes, unlike NGC300 X-l and 
IC10 X-l. The published PDS of M82 X-l and HoIX X-l 
do not extend to sufficiently low frequencies to determine 
whether they are also in the SPL state. We note that the 
broken power l aw P DS and power law PDS identified by 
ICropper et all (|2004h in NGC4559 X-7 and NGC4559 X-10 
are both artificial, caus ed by improper trea tment of non- 
sychronised lightcurves (jBarnard et alj|2007l ). 

5. Conclusions 

IC10 X-l is the first confirmed BH+WR system, and 
NGC300 X-l is thought to be one also. They are known 
to have strikingly similar X-ray periods and luminosities, 
so we studied them in detail. We have examined for the 
first time the PDS of NGC300 X-l, in four XMM-Newton 
observations, and also of IC10 X-l, in one XMM-Newton 
observation. We find that they are well described by a sim- 
ple power law with spectral index ~1, over the 0.002-0.1 
Hz range; such variability is characteristic of disc-accreting 
XBs in their high state, or of some HMXB pulsars that 
are powered by Bondi-Hoyle accretion. We also find a strik- 
ing resemblance between the Obs. 1 NGC300 X-l spectrum 
and the IC10 X-l spectrum. Since IC10 X-l is a confirmed 
WR+BH system, our results provide strong support for 
NGC300 X-l being one too. 

We have considered four accretion scenarios: disc accre- 
tion or Bondi-Hoyle accretion, onto a BH or NS. We favour 
disc accretion onto a BH for both systems, but cannot rule 
out Bond-Hoyle accretion onto a BH. The observed X-ray 
spectra of NGC300 X-l allowed us to rule out Bondi-Hoyle 
accretion onto a NS, but do not exclude disc accretion onto 
a NS+WR binary, or onto a NS LMXB that is merely co- 
incident with the WR. 

If both NGC300 X-l and IC10 X-l are indeed WR+BH 
binaries, then they are spectrally distinct from the known 
BH XBs. We propose that this difference may be due to 
these systems being persistently in the high state, allowing 
them to keep their disc corona, like LMC X-l; contrariwise, 
BH LMXBs enter the high state only in violent outbursts, 
where the corona may be ejected. NGC300 X-l and IC10 
X-l (and possibly LMC X-l) would then comprise a new 
class of BH binary, defined by their high mass transfer rates 
and consequent stable accretion disc configurations. This 
stable disc scenario may also explain why only ~10% of 
ULXs exist in the thermally-dominated high soft state. A 



fair comparison between the ULXs and NGC300 X-l and 
IC10 X-l requires detailed analyis of the variability of ULXs 
in addition to spectral analysis; we will carry this out in a 
future paper. 
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